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We present a method to design a finite decay rate for excited rotational states in polar molecules. 
The setup is based on a hybrid system of polar molecules with atoms driven into a Rydberg state. 
The atoms and molecules are coupled via the strong dipolar exchange interaction between two 
rotation levels of the polar molecule and two Rydberg states. Such a controllable decay rate opens 
the way to optically pump the hyperfine levels of polar molecules and it enables the application of 
conventional laser cooling techniques for cooling polar molecules into quantum degeneracy. 
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The fast decay of excited electronic state of atoms is a 
paradigm of interactions between light and matter. This 
decay is the main building block for several technolog- 
ical applications like laser cooling or optical pumping. 
In turn, these technologies open the way to the poten- 
tial implementation of quantum simulators and scalable 
quantum computers. In contrast to the atomic case, the 
decay of electronically excited states in polar molecules 
generally ends up in an unwanted excited rotational or 
vibrational internal state. This proves to be a major ob- 
stacle in combining the tools for controlling and cooling 
atoms with the strong electronic dipole-dipole interac- 
tions present in polar molecules. In this letter, we demon- 
strate the possibility to design an adjustable finite life 
time for an excited rotational level in a polar molecule, 
which overcomes the obstacles on the way towards optical 
pumping and laser cooling of polar molecules. 

Several alternative methods arc experimentally devel- 
oped in order to cool polar molecules which reside in the 
rotational and vibrational ground state to cold tempera- 
tures. For example, buffer gas cooling [1] and stark decel- 
eration [2]; see Ref. [3, 4] for a review. However, currently 
the most successful approach is to use ultra cold clouds 
of atoms from which molecules are created via photo as- 
sociation or coherent formation [5-7] . This technique has 
recently been perfectionized to create molecules in a well 
defined hyperfine state [8, 9]. Removing the remaining 
kinetic energy turns out to be extremely challenging due 
to the complex internal structure of the polar molecules 
[10]. Given these complications it may be profitable to 
use a hybrid setup of dipolar molecules with another sys- 
tem. A gas of Rydberg atoms is a promising candidate as 
they have both a large dipole moment and can be cooled 
to ultra cold temperatures [11-14]. 

In this letter, we demonstrate that a hybrid system 
of polar molecules with Rydberg atoms [15] allows for 
the design of a fast and controllable decay for an in- 
ternal rotational state. The main idea is based on the 
resonant and strong dipole-dipole interaction between 
two rotational states and two Rydberg levels: The polar 
molecules make a transition into the ground state, while 



at the same time the exchange of a virtual microwave 
photon excites the Rydberg atom into a higher level, see 
Fig. 1. Then, the spontaneous decay of the Rydberg level 
induces a finite effective decay rate for the rotational level 
of the polar molecules. 

We show that such a decay opens up the way for opti- 
cal pumping of hyperfine levels and laser cooling of polar 
molecules using standard techniques. Our setup has two 
main advantages. First, the decay rate is tunable to arbi- 
trary low values and, second, the decay takes place with 
very low recoil momentum. Consequently, the Doppler 
temperature is much lower than in their atomic counter 
parts and allows one to cool polar molecules into a regime 
where interactions and particle statistics play in impor- 
tant role. We demonstrate that the method exhibits the 
potential for laser cooling of strongly interacting many- 
body system into strongly correlated phases, such as the 
crystalline phase of polar molecules. We envisage, that 
our scheme is most suitable for polar molecules created 
by a coherent formation [8, 9], where unpaired atoms are 
naturally present, and the polar molecules are already at 
rather low temperatures. 
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FIG. 1: (a) Hybrid system: trapped polar molecules are in 
proximity to a cloud of a cold atomic gas with a separation 
I. (b) Relevant level structure: two rotational states for the 
polar molecule |e),; and \g}i are coupled via dipole-dipole in- 
teraction to two Rydberg levels \S)j and \P)j of an atom. In 
addition, the atom is driven from the ground state \G)j into a 
Rydberg level. The rotational level je) ; acquires a finite decay 
time due to the resonant coupling and the finite life time of 
Rydberg levels. 
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We start with the description of the setup for the en- 
gineering of a controUed decay rate for an internal rota- 
tional state of polar molecules. The decay takes place 
from the excited state denoted as je), into the ground 
state \g)i for the i"^ polar molecule, see Fig. 1. The tran- 
sition between the two states is characterized by a dipole 
matrix element d^g = (g|d|e) with d the dipole operator. 
The internal structure of the polar molecule is described 
by an effective spin 1/2 system with the Hamiltonian 
i?^^ = fii^o</2. The energy difference between the ex- 
cited and the ground state is TtWQ, and its typical value 
is in the range of the rotational splitting 2B. The polar 
molecule is coupled via the dipole-dipole interaction to 
the Rydberg levels of a close by atom. The relevant lev- 
els for the atom are the ground state \G)j and the two 
Rydberg states 15*)^ and \P)j- The ground state \G)j is 
coupled via a driving laser field to the Rydberg state \S)j. 
Within the rotating frame and using the rotating wave 
approximation, the Hamiltonian for the atom reduces to 



contains the polar molecules as well as the driven atomic 
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\G){S\,+\S){G\A, (1) 



where fl is the Rabi frequency of the driving laser field 
with the detuning h5 = Es — Eq. In addition, the Ryd- 
berg states are characterized by a dipole transition ma- 
trix element da ~ (5|d|P). 

The dipole-dipole interaction gives rise to a resonant 
exchange interaction between the polar molecules and the 
Rydberg levels \e)i\S)j — >■ \g)i\P)j- The setup is designed 
such that the exchange takes place near resonant with the 
detuning HA = Ep — Es — hajQ. The interaction between 
the atom and the molecule reduces to 
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where Vrf(r) = d^g ■ da/r^ - 3(deg • r)(da • r)/r^ denotes 
the dipole-dipole interaction with the positions of the 
polar molecules and Rj the coordinates of the atoms, re- 
spectively. The operators S~ = \g){e\i account for the 
transition within the polar molecule, while T^^ ~ \P){S\j 
is the corresponding operator acting on the Rydberg lev- 
els. 

The excited electronic states of the atom \S)j and \P)j 
are characterized by a finite life time with the decay rates 
7s and 7p, respectively. The time evolution of the full 
system is then described by the master equation for the 
density matrix p 

dtp = ^ [H, p] y [2Cap4 - clcaP - Pcica] , (3) 

a£{S,P} 

with the jump operators Ca = \G){a\ with a G {S,P} 
accounting for the decay from the Rydberg states into 
the ground state \G)j. In addition, the Hamiltonian H 



system, i.e., H = J2^ Hii> + H^J' + H,. 

The master equation Eq. (3) implies that any initial 
state ends up with all polar molecules in the rotational 
ground state \g)i, as can be seen by a numerical simu- 
lation of the time evolution. For weak couplings, it fol- 
lows that the dynamics is well described by an efi^ective 
decay rate for the excited rotational state of the polar 
molecules, which reduces to 
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(4) 



The last factor in the product ps = /{"fl/^ + 5"^) de- 
scribes the probability for an atom to be in the excited 
Rydberg state \S)j under the external drive . The second 
factor accounts for the finite line width of the Rydberg 
states and the first factor is the coupling matrix elements 
determined by the dipole-dipole interaction 
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Here, i/'m {4'jm) a-rc the initial spatial wave function and 
''Pin' {i'jm') the final spatial wave function for the polar 
molecule (atom), respectively. The decay rate in Eq. (4) 
is the main result of this paper. In the following we 
provide (i) a discussion of Eq. (4) and the limits of its 
applicability, (ii) numerical values relevant for concrete 
experiments, and (iii) a possible application in the from 
of Doppler cooling. 

Discussion of Eq. (4) — We start with analyzing 
the matrix elements My for two different setups: First, 
we consider molecules tightly confined in a three dimen- 
sional harmonic trap with oscillator length oq. If aq is 
much smaller than the distance between the atoms and 
the polar molecules Z, we can express the dipole-dipole in- 
teraction in terms of the harmonic oscillator wave func- 
tions. The excitation from the oscillator ground state 
is then suppressed by (ao/l)'^. Consequently, in lead- 
ing order the polar molecule remains within the spa- 
tial ground state. However, the resulting average energy 
transfer by 2h'^Tr'^ /mP corresponds to a recoil momentum 
of qd = 2'iTh/l. 

Alternatively, for polar molecules tightly confined in 
only two dimensions, see Fig. 1, a transfer of in-plane 
momentum K can appear. Its contribution is described 
by the in-plane Fourier transformation of the interaction 
potential Vd(Z,K) = / (ir|| Vrf(/, r|[)e'^''ii . It is peaked for 
K ^ l/l and decays exponentially for larger momenta. 
Consequently, also such a setup gives rise to a character- 
istic momentum recoil qd = 27rfi./L 

The above analysis reveals an important aspect of our 
setup: The recoil of the spontaneously emitted photon 
from the Rydberg state only acts on the atom, while 
the recoil for the polar molecule qd in turn is completely 
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determined by the separation I and can be strongly sup- 
pressed compared to the photon recoil. Note, that this 
discussion for the recoil energy is limited to the regime, 
where Fermi's Golden rule is applicable. Hence, we need 
to asses the validity of Eq. (4) . 

Let us return to the evaluation of Eq. (4). The de- 
cay rate is enhanced by the number of atoms within 
range of the dipole-dipole interaction, and can be writ- 
ten as J2j l-^iiP = The number of atoms 
within interaction range Ni strongly depends on the spa- 
tial distribution of the atomic sample and reduces to 
TV, = JdIi\Vd{n - R)p/9(R). Here, p(R) denotes the 
density of Rydberg atoms. For weak Rabi frequencies 
57/(5 <C 1 this density is determined by the atomic den- 
sity reduced by a factor ps accounting for the probability 
to find an atom in the Rydberg state. For strong drives 
the Rydberg-blockade phenomena starts to play a role 
and limits the Rydberg density; see [16] for a detailed 
discussion. In either case, this analysis provides us with 
an effective coupling strength 7Vi|Vd(l)p between a single 
molecule and the atomic cloud. 

We can now answer the question on the validity of 
Fermi's Golden rule which is based on weak coupling. In 
the presence of the enhancement by the different decay 
channels, this condition requires A//Vi'|Vd(l)| < 7p. As 
a consequence, the induced decay rate F^ is limited by 
the decay rate of the Rydberg level. However, the latter 
can be strongly enhanced and controlled by dressing the 
Rydberg level with a lower lying electronic state [15]. 
Denoting the decay rate of the additional electronic state 
by 7d, the far-detuned coupling of the Rydberg level \P)j 
with this electronic state with Rabi frequency fid and 
detuning provides a total decay rate for the dressed 
Rydberg level 7p ~ 7p + ^d/^^d- When providing 
numerical values for the obtainable decay rates for the 
rotational levels with take this enhancement into account. 

Numerical values — In the following, we provide ex- 
perimental parameters for LiCs polar molecules with the 
atomic reservoir given by Cs atoms excited into the Ryd- 
berg state. First, we ignore the weak hyper fine interac- 
tion for the polar molecules, and describe the character- 
istic energy scale; the incorporation of the nuclear spin is 
then straightforward. We choose for the excited states of 
the polar molecule the first excited rotational states with 
angular momentum one |e) = |1,0); here \j,m) denote 
angular momentum cigcnstates of the rotational degree 
of freedom. The dipole moment to the ground state re- 
duces to |deg| = d/^/^ with d — 5.5 Debye w 2.14eao the 
permanent electric dipole moment of the polar molecule. 
With a rotational splitting B « 5.84GHz, the transi- 
tion 1 1,0) \g) = |0,0) becomes nearly resonant for the 
Rydberg transition between states with principle quan- 
tum number n = 69 and (electronic) angular momentum 
I*) ~^ \P3/2)- According to the well established quan- 
tum defect theory [17] they have an energy difference of 
A « 40MHz and a dipole moment |da| « 10500 Debye. 



With these large dipole moments the resonant dipole- 
dipole interaction reduces to Vd{l) « 1.7 kHz even at a 
distance of / = 25/j,m between the polar molecules and 
the Rydberg atoms (separated along the z-direction) . 

This energy scale should be compared with the decay 
rate from the Rydberg state 7p, which is dominated by 
black body radiation and reduces to 7p ~ 840Hz. In the 
following, we will enhance this decay rate to 7p « 0.3MHz 
as discussed above. Consequently, we obtain a decay rate 
w 39Hz for each atom excited into the Rydberg level \S). 

Next, we need to estimate the enhancement Ni due to 
the coupling to a whole ensemble of Rydberg atoms. The 
characteristic blockade radius is in the range of « 3/xm 
[16], and consequently, the Rydberg density is limited by 
w lO^^cm"^; this density is much lower than the charac- 
teristic atomic density for a cold sample of Cs atoms. For 
an atomic sample within an oblate trap with transverse 
confining along the z-axis with 10/Ltm, the summation 
over the different Rydberg atoms provides an enhance- 
ment N, sa 52000. Then, the condition ^/NlVd{l) < 7p is 
satisfied and the application of Fermi's Golden rule for 
the estimation of the induced decay rate is valid. We fi- 
nally find a induced decay rate for the rotational level of 
the polar molecules F^ « 0.25MHz. 

Alternatively, it is experimentally interesting to work 
in a parameter regime with strong electric fields in com- 
bination with a strong confinement of the molecules into 
a two-dimensional setup [18]. However, the strong elec- 
tric field requires to work with much lower main principal 
number n of the Rydberg levels to prevent field ioniza- 
tion. 

In this setup the condition of resonant dipole-dipole 
interaction is conveniently obtained by selecting two Ry- 
dberg states which undergo a true level crossing within 
the stark map, but are coupled by exchanging a quan- 
tum of angular momentum. For Cs Rydberg atoms, 
the states with n ~ 17 and \S) = 1^5/2, ±3/2) and 
\P) = ](i5/2,±5/2) exhibit such a level crossing [19] at 
an electric field E w 3kV/cm. Moreover, they allow for a 
resonant dipole exchange with the two states \e) = Jl, ±) 
and \g) = 10,0) of the polar molecule [20]. On a dis- 
tance with I = 3/ini the dipole exchange energy reduces 
to Vd{l) « 3100Hz. With the decay rate for the Rydberg 
atoms "fs ~ Jp ~ 72kHz, the induced decay rate with a 
single Rydberg atom reduces to « 530Hz. Again, this 
decay rate is enhanced by the number independent de- 
cay channels Ni. Assuming the atomic system confined 
into a 2D layer above the polar molecules with a Ryd- 
berg density « lO^^cm^'^, the total decay rate reduces to 
r, « llkHz. 

Applications — An immediate application capitaliz- 
ing on the finite life-time of a rotational level in polar 
molecules is the optical pumping of the hyper fine struc- 
ture. The main mechanism is the electric quadrupole 
interaction, which couples nuclear spins to rotational ex- 
citations, and has been successfully used for the coherent 



4 



transfer of population between different hyperfine states 
[8, 9]. Note, that the optical pumping may give rise to co- 
herence between polar molecules separated on distances 
smaller than I; it is this coherence, which allows to apply 
conventional laser cooling techniques [21-23] to temper- 
atures below quantum degeneracy. 

In the following, we describe the peculiarities of our 
setup for Doppler cooling. However, we would like to 
stress, that the induced decay rate can also be used for 
more sophisticated cooling techniques, such as Sisyphus 
cooling, which in general will give rise to more efficient 
cooling of the system. In order to achieve Doppler cool- 
ing, the transition into the excited rotational state is 
driven via a Raman transition with Rabi frequency fl, 
the detuning from resonance A, and a momentum trans- 
fer q. In contrast to conventional laser cooling, here, 
the spontaneous decay from the excited rotational state 
transfers a different momentum to the polar molecule. 
Consequently, we have to distinguish between two mo- 
menta. The momentum q absorbed in the Raman pro- 
cess and the small momentum qd from the spontaneous 
decay; see above for a discussion of this momentum. 

The standard laser cooling techniques are character- 
ized by two energy scales: the higher energy /temperature 
= mr^/(2q^) determines the kinetic energy, where 
laser cooling becomes effective; here F denotes the in- 
duced decay rate from the excited rotational state. Sec- 
ond, the equilibrium state induced by the laser cooling is 
characterized by the Doppler temperature Td = hT. 

A major advantage of the present scheme is the pos- 
sibility to tune the recoil momentum q independently 
of the decay rate F. As a consequence, the capturing 
temperature and the Doppler temperature can be ad- 
justed for optimal cooling. The Doppler temperature of 
the equilibrium state can be continuously decreased by 
reducing the decay rate F. The lowest temperature is 
therefore only limited by the small recoil momentum 
from the spontaneous decay, i.e., T^^^ = h'^q'^/2m. This 
temperature is in general much lower than the Doppler 
temperature for conventional laser cooling of alkali atoms 
[23], and the present scheme allows to laser cool polar 
molecules to temperatures comparable to quantum de- 
generacy. 

Note, that for conventional Doppler cooling, the re- 
absorption of emitted photon provides a decrease in cool- 
ing efficiency for dense atomic samples. The analogue ef- 
fect in the present scheme would be the re-excitation of a 
polar molecule by a resonant dipolar exchange. However, 
this rate is strongly suppressed for a fast decay of the Ry- 
dberg level and a proper balance between the number of 
atoms N and polar molecules M with the requirement 
M < TV. 

Close to quantum degeneracy, the statistics of the par- 
ticles starts to play an important role. The Doppler cool- 
ing scheme described above remains valid also in this 
regime. This allows one in principal to cool directly into 



a Fermi sea for fermions or a Bose-Einstein condensate 
for bosons. However, the Doppler cooling scheme breaks 
down at low temperatures and high particle densities, 
where interactions between the molecules become rele- 
vant. 

There are two main effects due to interactions: first, 
the interaction potential between two molecules in state 
\g), and the interaction potential between the states \g) 
and |e) are in general different; especially, as for the two 
states \g) and |e) also a resonant exchange of a microwave 
photon can appear. During the Doppler cooling pro- 
cess, each quantum jump with the excitation of a po- 
lar molecule into the higher rotational state provides a 
fluctuating force to the surrounding polar molecules and 
eventually provides a source of heating. Second, the dif- 
ference of the interaction potentials leads to a level broad- 
ening of the excited rotational state, and the desired en- 
ergy selectivity for the Doppler cooling breaks down as 
soon as the level broadening reaches the induced decay 
rate. 

However, one can envisage a strategy that overcomes 
this limitation and potentially allows cooling directly into 
a strongly interacting many-body phase such as the crys- 
talline state [18]. The trick is to adiabatically eliminate 
the excited rotational states by irradiating a far-detuned 
Raman laser. Consequently, the many-body states ac- 
quire a finite lifetime inherited from the small admix- 
ture of the excited state. In the following, we denote 
by I a) and |/3) eigenstates of the fully interacting many- 
body system with energies Ea and Ep, respectively. Us- 
ing Fermi's Golden rule, we obtain transition rates Ta^p 
from the initial state \a) into the final state \/3) due to 
the coupling of the system to the bath 
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Here, F denotes the decay rate of the bath, while the 
matrix elements for the many-body system reduces to 



Map ^ / dx (a|i/''i'(x)'0(x)e' 



jqxi 



(7) 



It is important to point out, that the energy dependence 
is in analogy to conventional laser cooling techniques, i.e., 
for Ea — Ep > and blue detuning A > 0, the transitions 
\a) —T' 1/3) take place with a higher rate than the reversed 
transitions — > \a), and consequently leads to a cooling 
of the system. 

Conclusion. In summary, we proposed a scheme to 
induce a controlled decay rate for rotational levels of po- 
lar molecules. We showed how one can harness this finite 
lifetime to perform Doppler cooling into quantum degen- 
eracy. Finally, we commented on how to implement a di- 
rect cooling into a strongly interacting many-body phase. 
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